Enhanced stability of tetratic phase due to clustering 
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We show that the relative stability of the nematic tetratic phase with respect to the usual uni- 
axial nematic phase can be greatly enhanced by clustering effects. Two-dimensional rectangles of 
aspect ratio n interacting via hard interactions are considered, and the stability of the two nematic 
phases (uniaxial and tetratic) is examined using an extended scaled-particle theory applied to a 
polydispersed fluid mixture of n species. Here the i-th species is associated with clusters of i rect- 
angles, with clusters defined as stacks of rectangles containing approximately parallel rectangles, 
with frozen internal degrees of freedom. The theory assumes an exponential cluster size distribu- 
tion (an assumption fully supported by Monte Carlo simulations and by a simple chemical-reaction 
model), with fixed value of the second moment. The corresponding area distribution presents a 
shoulder, and sometimes even a well-defined peak, at cluster sizes approximately corresponding 
to square shape (i.e. i ~ «), meaning that square clusters have a dominant contribution to the 
free energy of the hard-rectangle fluid. The theory predicts an enhanced region of stability of the 
tetratic phase with respect to the standard scaled-particle theory, much closer to simulation and to 
experimental results, demonstrating the importance of clustering in this fluid. 

PACS numbers: 61.30. Cz, 61.30.Hn, 61.20. Gy 



I. INTRODUCTION 

The hard rectangle (HR) fluid constitutes a paradig- 
matic example of a two-dimensional fluid exhibiting 
surprisingly complex phase behavior: different phase 
symmetries, phase transitions with different order, and 
defect-mediated continuous transitions of the Kosterlitz- 
Thouless type all governed solely by entropy. This 
peculiar two-dimensional system has three equilibrium 
fluid phases: isotropic (I), where particle axes are ran- 
domly oriented, uniaxial nematic (N u ), with particles 
preferentially aligned along a single nematic director, and 
tetratic nematic (Nt), possessing two equivalent perpen- 
dicular nematic directors, with long particle axes oriented 
along one of two directors with equal probability. 

In a pioneering study, Schlacken et al. Q applied 
scaled-particle theory (SPT) on a fluid of HRs to demon- 
strate the stability of the N t phase, a phase which cannot 
be stabilised in a fluid of hard ellipses. The intersection 
between the two spinodals associated with the I-Nt and 
I-N u transitions defines a limiting aspect ratio n = L/ao 
(with L and oq the length and width of the rectangles, 
respectively) for the stability of the Nt phase, which is 
located at n ~ 2.62. Thus, for lower values of aspect 
ratio, the isotropic fluid exhibits a continuous transition 
to the Nt phase, whereas if k > 2.62 the I phase goes 
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directly to the uniaxial nematic phase N u . 

The study of Schlacken et al. was later supplemented 
by the calculation of the complete phase diagram [3], also 
within the context of SPT. It was found, in particular, 
that the Nt phase undergoes a transition to the N u phase 
at high density, the nature of which changes from second 
to first order at a tricritical point. In addition, it was 
shown 01 that, at the level of a particular approximation 
for density-functional theory, the Nt fluid is metastable 
with respect to a phase with (either partial or complete) 
spatial order; the theory was approximate in the sense 
that it included the exact functional form of two-body 
correlations, but only approximate higher-order correla- 
tions. On the other hand, Monte Carlo (MC) simula- 
tions conducted on hard squares [4j and on a HR system 
[B[ with k = 2 indicated, as expected, that the fluid ex- 
hibits quasi-long-range tetratic order, and that the high- 
density phase consists of an aperiodic crystalline tetratic 
phase exhibiting random tiling on a square lattice. MC 
simulation of a HR fluid confined in a slit pore Q show 
the presence of weak tetratic correlations at the centre of 
the pore. 

In the experimental front, recent results for colloidal 
discs forced to stand on edge by external potentials [l[ 
(and hence interacting approximately as HRs) have also 
demonstrated that tetratic correlations play a vital role 
in this system. Also, experiments conducted on a mono- 
layer of vibrated granular cylinders lying on a plate Q 
have shown tetratic correlations for cylinders with aspect 
ratio as high as n = 12.6. 

In Ref. strong evidence, based on MC simulations, 



was presented for the thermodynamic stability of a Nt 
fluid when k is at least as large as 7. These simulations 
were supplemented by an extended SPT model that ex- 
actly incorporates the second and third virial coefficients 
while resumming the remainder of the virial series. The 
inclusion of the third virial coefficient increases the in- 
terval in k where the N t phase is stable, approaching the 
simulation result. Specifically, the I-Nt transition line 
moves to lower packing fractions, while the intersection 
point between the I-N u and the I-N t transitions shifts to 
k = 3.23. However, this value is still lower than the value 
indicated by simulations. 

The properties of the HR fluid are to be contrasted 
with those of hard discorectangles. MC simulations of 
this system have been conducted [§], and the global phase 
diagram was computed. A careful inspection of particle 
configurations in the isotropic phase shows that there are 
peculiar equilibrium textures, with large clusters contain- 
ing particles arranged side by side, exactly as in the HR 
fluid. These configurations are favoured by the particular 
shape of the particles and by the reduced dimensional- 
ity. However, in contrast with the HR fluid, neighbouring 
clusters do not exhibit strong tetratic correlations and, 
therefore, the formation of a tetratic nematic phase is 
discouraged in the hard discorectangle fluid. 

In this article we address the problem of how the 
present theoretical understanding of the HR fluid can 
be improved by consideration of clustering effects. Our 
thesis is that these effects, very apparent in our own MC 
simulations but not addressed by the theories proposed 
up to now, are a key factor in the stabilisation of the 
Nt phase. Inclusion of cluster formation is responsible 
for the enhancement of the region of N t stability in the 
phase diagram. In the model proposed, clustering is ap- 
proximately taken care of by treating clusters as distinct 
species in a mixture of polydispersed rectangles. The 
functional form for the cluster size distribution is as- 
sumed to be exponential, an assumption supported by 
cluster statistical results based on MC simulations and 
by a simple chemical-reaction model (see Appendix) , and 
is introduced in the model as an input. The thermody- 
namics of the polydispersed mixture is analysed using 
SPT. The results indicate that clustering (assimilated in 
the theory by means of a polydispersity parameter) sta- 
bilises the N t phase for values of aspect ratio much higher 
than k — 2.62 if the polydispersity is sufficiently high. 
Polydispersity parameters obtained from simulation give 
support to the model. 

The article is organized as follows. In Section II 
we present numerical evidence that the size distribu- 
tion in the HR fluid is an exponentially decaying func- 
tion. Section III presents the main ideas of the theo- 
retical model proposed. The conclusions are drawn in 
Section IV. Finally, details on the model and on the pro- 
cedure of solution are relegated to Appendices A and B, 
while Appendix C contains a chemical-reaction model for 
monomer aggregation which also supports the assump- 
tion of an exponential cluster size distribution. 
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FIG. 1: (Colour online). Configuration of hard rectangles 
of aspect ratio K — 7 at packing fraction r\ — 0.655, as ob- 
tained from MC simulation. Clusters, defined by the pair con- 
nectedness criterion explained in the text, have been coloured 
according to their size. 



II. MONTE CARLO SIMULATION OF 
CLUSTERING 

We started by applying standard isobaric (NPT) MC 
techniques on a two-dimensional fluid of hard rectangles, 
with aspect ratios n — 3, 5 and 7, using N = 1400 rect- 
angles. The transition from the I phase to the Nt phase 
was identified approximately by inspection of the tetratic 
order parameter 92 = (cos 40), where is the angle be- 
tween the long axis of the particle and an axis fixed in 
space. Once the samples were equilibrated, cluster statis- 
tics was applied. The criterion for pair connectedness, 
i.e. for deciding when two neighbouring rectangles can 
be considered to be 'bonded', was based on the relative 
angle 4>n between the long axes of the particles and their 
relative centre-of-mass distance ri2, by demanding that 
012 < 5 and | r 12 1 < e. Typical values adopted were 
S = 10° and e = 1.3tro, although the conclusions to be 
presented below do not seem to depend qualitatively on 
the exact values (provided they are not too large). 

Fig. [TJpresents a configuration of rectangles with n = 7 
at a packing fraction rj = 0.655. This corresponds to a 
tetratic phase. Identified clusters have been coloured ac- 
cording to their size. One can see how these clusters look 
like large "super-rectangles" arranged along two perpen- 
dicular directions. Therefore, the tetratic structure is 
maintained not only for single rectangles, but also at 
the level of clusters ("polydispersed super-rectangles"). 
This hierarchical feature of the tetratic symmetry will 
give support to the theoretical model to be presented in 
the following section. In Fig. [2] a logarithmic histogram 
of the size distribution (averaged over configurations) is 
shown; in the figure, Xi is defined as the fraction of clus- 
ters of size i (see next section). The distribution looks 
exponential. All the curves in the figure pertain to a 
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FIG. 2: Size distribution function obtained from MC 

simulation, in logarithmic scale, for a fluid of HRs with differ- 
ent values of aspect ratio and packing fraction: filled circles, 
k = 3 and r\ = 0.635; open squares, k = 7 and r\ = 0.600; and 
grey squares, k = 5 and r\ = 0.657. Straight lines are linear 
fits. Inset: area distribution function ixi for the case k = 5 
and r) = 0.692. 
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FIG. 3: Monomer h m (4>) (solid line) and cluster h c (<f>) 
(dashed line) orientational distribution functions for the case 
k = 7 and r\ = 0.655, as obtained from MC simulation. 



property does not seem to be followed in the uniaxial 
nematic phase. 



tetratic phase, but the same behaviour is observed also 
in the isotropic phase (with the general trend that, for 
given ac, the slope decreases, i.e. the size distribution 
function becomes wider, hence the average cluster size, 
as density increases). In the size region corresponding 
to square clusters (i.e. clusters with aggregation num- 
ber i ~ k) it is possible to see an incipient shoulder that 
grows as density is increased. This feature is more appar- 
ent in the area distribution function, ixi (giving the frac- 
tion of area occupied by clusters of a given size), which 
presents a shoulder and sometimes even a well-defined 
peak (inset of Fig. [2]), indicating that square clusters are 
structurally very relevant and contribute very decisively 
to the thermodynamic properties of the fluid. 

We end this section by discussing the orientational dis- 
tribution functions. Let h m ((j)) be the monomer orien- 
tational distribution function, giving the probability of 
finding a given rectangle with its long axis forming an 
angle with respect to the director. Having defined clus- 
ters in the fluid, we may also define an orientational dis- 
tribution function associated with clusters, h c (</>), giving 
the probability of finding an average cluster (regardless 
of its size) oriented with angle 4>. In the simulations we 
compute h c (4>) by averaging over all the identified clus- 
ters, and then over all MC configurations. Fig. [3] gives 
the functions h m (<fi) and h c {4>) for a state with tetratic 
symmetry. It is interesting to note that these functions 
are almost identical (with the same symmetry and, con- 
sequently, with peaks with the same height), so that 
the structure of tetratic ordering is maintained from the 
monomer level to the cluster level, the functions h m (4>) 
and h c (4>) obeying a kind of "similarity" property. This 



III. THEORETICAL MODEL 

The model we propose accounts for clustering effects in 
an approximate way. The clear identification of clusters 
in the simulations, along with their approximate rect- 
angular shape, leads to a simple model where clusters 
are regarded as single rectangular particles with no in- 
ternal degrees of freedom. Therefore, we consider a n- 
component mixture of two-dimensional hard rectangles 
of dimensions L and cr, , with Oj = i<jQ and i the aggrega- 
tion number. Each of these rectangles is assumed to be 
composed of i monomers in perfect contact in a side-by- 
side configuration. The total number of monomers, Nq, 
can be written as 

n 

N = J2i^h (1) 

where Ni is the number of clusters containing i 
monomers. Dividing by the volume V , 

n 

p = X/^> ( 2 ) 

where po = Nq/V is the total monomer density, and 
Pi = Ni/V = Xip is the density of clusters of size i, with 
Xi = Ni/N their number fraction, while N = J^i ls 
the total number of clusters and p — N/V their den- 
sity. The set {x{\, i = 1, • • • , n, is a central quantity in 
our model, since it contains information about clustering 
tendencies. We will assume Xi to be an exponential with 
i, as explained later. The total packing fraction of the 
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FIG. 4: (Colour online). Multicomponent isotropic mix- 
ture composed of different species, where species, depicted in 
different colours, correspond to clusters of a particular size. 
Clusters are defined as aggregates of rectangular monomers 
in a side-by-side configuration. In this instance the monomer 
aspect ratio is k = 5. 



system is 



(3) 



where dj = ia (a = Lao) is the area of a cluster of size 
i. Using ([2]), we easily obtain r\ — poaQ. 

Fig. |4] shows a schematic representation of an isotropic 
configuration of clusters of different sizes (with monomer 
aspect ratio n = 5) . The original monomers that give rise 
to each cluster are indicated but note that, in our model, 
the identity of the monomers is lost, as monomers in the 
same cluster do not interact dynamically, always being 
in perfect side- by-side contact. 

The free-energy of this multicomponent mixture of 
hard rectangles will be modelled by means of SPT ap- 
plied to a fluid mixture of freely-rotating hard rectangles. 
The orientational properties of the mixtures will be char- 
acterised by orientational distribution functions hi (4>) for 
each component i. The free-energy density functional 
<f> = (3F/V is written as 



*KM] = p E : 



ln^Vj) 



- l + lny + yS[{hi}} 




dcj>hi{4>) ln[hi((j>)ir] 



(4) 



where Vj is the thermal volume of i-sized clusters, and we 
defined y = p<io/(l — v)- Due to the head-tail symmetry 
of the particles, the angle <f> can be restricted to the inter- 
val [0, 7r] and the functions hi(<f>) normalised accordingly. 



The function S [{hi}] 



J2t j XiXjSij [{hi}], with 



Sij [{hi}] 



(n + ijn ((I sin </) tj \)) 



+ 5^ + Ml cos ^-1)), 



(5) 



is related to Aij, the angle-averaged excluded area 
between clusters i and j, as Sij — (Aij /oq — i — 
2)1% The shorthand notation ((f(<ftij)}} has been used 
for the double angular average of a generic function: 
((f(4>ij))) = Jo d4i d(/)jhi((l)i)hj((t>j)f((j)ij). Now a bi- 
furcation analysis of (dJ at the I-N U) t transition (see Ap- 
pendix) allows us to obtain the packing fractions of the 
TN U] t spinodal lines as 
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with k 



1 for the uniaxial and k 

,(«) 



2 for the tetratic ne- 



matic, while m "' — ^ Xii a (a — 1,2), are the first and 
second moments of the cluster size distribution function. 

Based on the MC results, we adopt an exponential clus- 
ter size distribution: 



1-g 
l-q n 



1,- 



(7) 



with q = e~ A (A > 0). The prefactor in ([7]) ensures that 
the distribution is normalised, i.e. that Yli x i = !• The 
first two moments can be derived analytically: 



(8) 
(9) 
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(l + n(l-g)) 2 
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q?(l-q n ) 





Now we use §§§ to obtain the maximum aspect ratio 
which can support a stable tetratic phase. This follows 
by imposing the condition rj^ t = r/g u , i.e. by searching 
for the intersection point of the two spinodal lines I-N u 
and I-Nt in the phase diagram r\ — n. Solving for the 
corresponding value of k, we obtain 



3m { y ± a/9 



7/7 



(2) 



(10) 



In the specific case where the number of species goes to 
infinity, n — > oo, we obtain from f8l)- lfT0|) 



- (2k 2 - 3« - 1 ± \fn 2 + 6k + 1 



-1 1/2 



,(H) 



where A = \ /wq 2 '/ ^m^'j — 1 = y/q, a measure of 

polydispersity, is the relative mean square deviation. The 
two functions A(k), corresponding to the two signs in 
(fTTjl . are plotted in Fig. [5] (note that one of the branches, 
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FIG. 5: The functions A(/t) in terms of inverse aspect ra- 
tio Symbols indicate the nematic phase which is stable 
in the region in question. The region of tetratic stability is 
shaded. Open squares indicate values of polydispersity, for 
values of aspect ratio n = 3, 5 and 7, at the I-Nt transition, 
as obtained from MC simulation; error bars correspond to 
uncertainty in A originating from uncertainty in location of 
spinodal. 



the one with minus sign, can be only calculated for 
k > n* = (3 + \/5)/2). The two lines define a region 
(shaded in the figure) where the tetratic phase may be 
stabilised. Thus we see how, as the value of polydisper- 
sity A is increased, the maximum value of k for which 
the tetratic phase ceases to exist increases, which means 
that polydispersity enhances the formation of tetratic or- 
dering. 

Note that the values of A at the upper boundary of 
the tetratic region are very high, which means that there 
should be a large number of clusters with their long side 
ia » L, while their width is L. This causes the effective 
average cluster aspect ratio iao/L to be very high, which 
induces formation of uniaxial nematic order and thus the 
Nt phase is destabilised in favour of the N u phase. The 
presence of such big clusters is not observed in the sim- 
ulations, which means that a more sophisticated model 
should somehow include additional entropy terms that 
discourage the formation of big clusters; for the sake of 
simplicity, we have kept the ingredients of our model to a 
minimum and avoided any such additional complications. 

Finally, in Fig. Rvalues of polydispersity A calculated 
from the simulated cluster size distribution function have 
been indicated by symbols. These values correspond to 
the estimated spinodal line of the I-N t transition for the 
cases k — 3, 5 and 7. Since these estimations are very 
rough, the uncertainty in packing fraction at the spin- 
odals goes over to the value of polydispersity (which nat- 
urally depends on k and rj), which is represented in each 
case by error bars. We can see that in two of the cases 
the symbols are well inside the tetratic stability region 



calculated from the theory. In the case n — 7, where 
uncertainty is larger, the MC estimate of A lies outside 
(but close to) this region. 



IV. CONCLUSIONS 

Because of the reduced dimensionality, fluids of two- 
dimensional hard anisotropic particles exhibit strong 
clustering effects: particles have less freedom to orient in 
space, which fosters configurations where neighbouring 
particles lie parallel to each other. However, the impact 
of this on the onset of new macroscopic symmetries de- 
pends very sensitively on the particular geometry of the 
particles. Thus, in fluids of rectangles, neighbour clus- 
ters have a strong tendency to adopt orthogonal relative 
configurations, since these clusters are almost perfect big 
rectangles made of several, almost parallel, monomers. 
These strong tetratic correlations are capable of generat- 
ing full macroscopic tetratic order and a thermodynam- 
ically stable tetratic phase. Therefore, clustering effects 
are crucial to understand phase behaviour in the HR fluid 
(and possibly also in the HDR and related fluids), but 
simple theories at the level of two-body monomer cor- 
relations (Onsager, SPT, etc.) cannot account for these 
effects. 

In this work we have presented a simple theory that 
incorporates clustering in terms of cluster polydispersity, 
where clusters are considered to be inert particles with no 
internal degrees of freedom. This assumption may be ac- 
curate provided the cluster lifetime (i.e. the average time 
it takes for a cluster to disappear since it was formed) is 
longer than the typical cluster diffusion rates in the fluid. 
Validation of this condition will have to wait for molec- 
ular dynamics simulations of the hard-rectangle fluid. 
Once the fluid is modelled in terms of a multicompo- 
nent mixture, one of the available theories for mesophase 
formation can be used. We have used SPT and have ex- 
amined the consequence of polydispersity in the phase 
diagram. As expected, polydispersity enhances the sta- 
bility of the tetratic phase. Due to the limitations of our 
model (e.g. the cluster size distribution has to be im- 
posed from outside and does not result from the theory) , 
we cannot make any quantitative comparison with avail- 
able simulation and experimental results. However, the 
model can qualitatively explain the formation of tetratic 
order for rather high values of aspect ratio, as shown by 
simulation and experiment. 



Acknowledgments 

Y.M.-R. gratefully acknowledges financial support 
from Ministerio de Educacion y Ciencia (Spain) under 
a Ramon y Cajal research contract and the MOSAICO 
grant. This work has been partly financed by grants 
Nos. FIS2005-05243-C02-01 and FIS2007-65869-C03-01, 
also from Ministerio de Educacion y Ciencia, and S- 



6 



0505/ESP-0299 from Comunidad Autonoma de Madrid 
(Spain). 



non-linear equations which, after some algebraic manip- 
ulations, can be transformed into a set of equations for 
the new variables 



V. APPENDIX 

In this appendix we provide additional details and fur- 
ther information on the consequences of the model. It 
contains two sections. In Section A, details on the set 
of non-linear equations that have to be solved to obtain 
the equilibrium properties of the HR fluid are provided. 
Also, the bifurcation analysis of the I-N u ,t transitions is 
presented. Section B is devoted to discussing the nature 
of the different phase transitions, together with the be- 
haviour of the distribution functions and to a comparison 
with simulations. In Section C, a simple chemical model 
of aggregation is discussed. 



A. Minimisation of free energy and bifurcation 
analysis 

Using Fourier series to represent the orientational dis- 
tribution functions, 



hi((f>) = -J2 h k ] cos(2fc0), 



(12) 



fc>0 



with = 1 Vi, together with Eqn. ((5]), we find 



1J 



*+-r + (-ir(<+j) 



rf4 J,) > (13) 



where g k = -(1 + S k0 )/2(4k 2 - 1). Defining 
rn^ =Y,*iV*hf , a = 0,1, 

i 

we obtain 

^XiXjSij = ~^2gkS 2 k , s k = m k 0) + (-1) 



(14) 



(i) 
k m k 



(15) 



Note that m, 



(0) 



J2i X i = 1 While = J2i x i i is 

the first moment of the discrete cluster size distribution 
function {xi}. Using this notation, the free-energy per 
particle tp = §/po can be written 



LP = In — - 1 + ^2 X l\ m («iVt) 



i=l 



dcj>hi{4>) In [hi(<j))ir] 



(i) 

nm" k 



TjJ2 gkS k> ( 16 ) 



with yo — 77/(1 — 77). The functional minimization of 
(fT6|) with respect to hi{4>) gives a set of self-consistent 



s k 



'■k ' 



qW = / d4>cos(2k(t>)h i ((f>), 



(17) 
(18) 







where the normalized orientational distribution functions 

are 



hi((/>) 



4j/ok \ - 

— w 2- Sfc 

k>l 



(19) 



(-i) fc ; 



g k cos(2/c0). 



(20) 



The linearization of l|17p with respect to Sk (k > 1) allow 
us to obtain the expression © for the packing fractions 
at the I \ u .f spinodal lines. 



B. Phase transitions and distribution functions 

In this subsection we analyse the free energy branches 
of the model in order to understand the nature of the 
different phase transitions. The cluster distribution func- 
tion Xi of the mixture is assumed to be exponential, and 
the width of the distribution is fixed via the polydisper- 
sity parameter q (or A). Eqns. fT7j) - (|20"ll are solved 
for different values of 77 to find all metastable and sta- 
ble phases, either I, N u or N t phases. In Fig. [6] the 
free-energy branches as a function of 77" 1 for different 
values of q (and hence for different polydispersities) are 
shown. As can be seen, the Nt phase begins to be sta- 
ble from q w 0.30. Also, it is clear that, for q — 0.25 
and 0.35, the I-N u or N t -N u transitions are of first order 
(free-energy branches cross with different slopes). The 
coexistence values of 77 cannot be determined from the 
standard double-tangent construction, since the present 
system is polydisperse. The usual procedure then is to 
fix the distribution function xf 1 and packing fraction 77'°' 
for the parent phase (I or N U)t phases) and find the cloud 
and shadow curves. We have not implemented this pro- 
cedure here. However, for those values of q for which the 
transitions are continuous (i.e. q — 0.5 and 0.65), the 
present procedure adequately determines the transition 
densities. A similar situation occurs for the cases n = 5 
and k = 7 (not shown). 

It is also interesting to look at the orientational dis- 
tribution functions of monomers and clusters. From the 
corresponding functions for clusters of size i, i.e. /ij (</>), 
it is easy to define a cluster orientational distribution 



(c) 



(d) 



FIG. 6: Free energy per particle <3?/p vs. inverse packing 
fraction 77 _1 from SPT results for the multicomponent HR 
fluid of monomer aspect ratio n = 3 and with different val- 
ues of polydispersity: (a) q = 0.25, (b) 0.35, (c) 0.50 and 
(d) 0.65. Continuous curves: tetratic phase; dashed curves: 
uniaxial nematic phase; dotted curves: isotropic phase. In (c) 
and (d) symbols indicate bifurcation points from the isotropic 
to the tetratic (filled circles) and uniaxial nematic (open cir- 
cles) phases. In (a) and (b) a straight line in rf~ 1 has been 
subtracted to better visualise the curves. 



function h c (4>) as 



(21) 



i=l 



From this, order parameters of the multicomponent mix- 
ture can also be defined: 



Q 



(fe) 



i(fc) 



k = 1,2. 



(22) 



In the case of monomers the situation is a bit more 
complicated, since in our model we have lost track of 
monomers as distinct entities. However, we can simply 
count the number of monomers pointing along some an- 
gle <p from the set hi(<j>) and then divide by the average 
number of clusters. Here we have to bear in mind that 
clusters of size i (having i monomers) with i < n (and k 
an integer) have a long axis in a direction perpendicular 
to that of clusters with i > k: therefore we write: 



h m {4>) = 



1 



id 



(i) 



}^Xjihj((f)) + 2J Xiihi(<j) + ir/2) 



i=n + l 



(23) 



FIG. 7: Cluster (a) and monomer (b) orientational distribu- 
tion functions for k = 5 and q = 0.65. The continuous line 
corresponds to a stable N u phase with r\ = 0.85, whereas the 
discontinuous line is for the metastable Nt phase at the same 
value of r\. 



parameters of the monomers can be calculated as 



<cos(20)) r 



(cos(4(/>)) 7 



E ■ 

i=no + l 



j no 

— r 2J XiiQi ■ 



(24) 
(25) 



'0 i=l 



Even for values of aspect ratio for which there exists a re- 
gion of tetratic stability, the N u is always the more stable 
phase for high values of packing fraction. It is interesting 
to note that, in this situation, the cluster and monomer 
distribution functions h Cim (4>) usually have a secondary 
peak at (f> = 7r/2 corresponding to tetratic ordering (this 
feature is also present in the simple SPT for the one- 
component HR fluid). Fig. [7| shows these distributions 
for the case k — 5 and q = 0.65. Note that, according to 
Fig. \E[ there exists a stable tetratic phase in this case, 
with a value of packing fraction at the N t -N u transition 
of 77* = 0.839. In the figure, the packing fraction cho- 
sen is 77 = 0.85 > rf and, therefore, the stable phase 
has an orientational distribution function pertaining to a 
N u . However, both the cluster and the monomer orienta- 
tional distribution functions of N u have secondary peaks. 
We have also plotted in Fig. [7] the distribution functions 
of a metastable Nt at the same value of 77. Finally, it 
is also interesting that, although the distribution func- 
tions h m {<j>) and h c (4>) in the tetratic region are always 
tetratic-like (i.e. all maxima have the same height, as it 
should be by construction), they do not coincide in the 
uniaxial nematic phase and, consequently, the similarity 
property at work in the tetratic phase is not obeyed for 
the uniaxial nematic (see Section [TT}l . 



C. Chemical reaction model 



with n = [k] (note that in the case where k is not an 
integer this division has to be done also). Thus the order 



Aggregation phenomena in dilute fluids (e.g. micelle 
aggregation) are very often described in terms of a chem- 
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Eqns. l(26|) . together with ([28| . give ^ = xiq 2 1 with 



80 100 



FIG. 8: Packing fractions r\ of the I-Nt (continuous line) and 
I-N u (dashed line) transitions versus aspect ratio n. 



ical reaction model. An exponentially decaying size dis- 
tribution immediately emerges from these models. Here 
we exploit the idea and carry it further using our density- 
functional approximation. This model assumes that the 
lifetime of a cluster is sufficiently long that it can be de- 
fined as a distinct 'chemical' species. 

One assumes a chemical reaction of the type C/ + Ci # 
C/ + i, with Ci denoting a 'chemical' species (i.e. a cluster) 
containing I monomers. Chemical equilibrium between 
clusters and monomers then implies the relation 



m = ifJ.i, i = 2, • • • ,n. 



(26) 



Now the chemical potential of the i th species can be 
calculated from our density-functional theory as 



01* = (27) 
dpi 

which results in 

= In (xjVi) + / d(j)hi(4>) MM0) 7r ] + toy 
Jo 

+yi + 2y^2 x j S ij + V 2si - (28) 

3 

The n—1 Eqns. lj26|) . together with the condition (jJJ, are 
a set of n equations with n unknowns (x2 , • • • , x n and p) 
which allow to find the equilibrium configuration of the 
fluid. Due to the simplicity of the model an analytical 
solution can be found. 

Since, for the isotropic phase, we have 



Sij) = 



(29) 



yo 



2yo 



+ 1 



(30) 



Here we assumed that V, = V{ (consistent with the ab- 
sence of internal degrees of freedom in the clusters and 
also with the assumption that all clusters have the same 
mass). Now, since J2i Xi — 1, we find, for a fluid with an 
infinite number of species, 1 = Xi/(l—q) which, together 
with (30]), give 



xi = <l 



yo 

— exp 




(31) 



Also, the first moment can be calculated self-consistently 



as 




(32) 



This solution means that Xi = (l — l/m^^J /m^K 

While the first moment is the solution of Eqn. (J32j) , 
the second moment results in 



,(2) 



lEl2jr<; 

i=l 



(i- q y 



to 1) (2to 1) - 1). 



(33) 



(i) 



The polydispersity coefficient, defined as A 

Now the I-N u and I-Nt spinodals can be calculated by 
solving Eqn. ([32]). together with the value of yo obtained 
from Eq. Jfj]). Fig. [8] contains the functions rj u ,t{^) ob- 
tained as the solutions of f32|) and ((6]). As can be seen the 
uniaxial nematic is more stable than the tetratic up to 
k* sa 41.65. This results from the peculiar behaviour of 
the area distribution function ix%. its zeroth-order mo- 
ment is always greater than k, while it is larger in the 
tetratic phase when k < k* [Fig. [U(a)]; this behaviour 
is inverted for k > k* , and the moment becomes larger 
for the uniaxial nematic phase, as can be seen in Fig. 
[9|b). Again this peculiar behaviour is due to the rela- 
tively high proportion of very big clusters, with iao > L, 
which stabilise the N u phase against the N t phase. This 
behaviour is not observed in simulations, because the for- 
mation of very big clusters is penalised by fluctuations. 
The model has the value that an exponential cluster dis- 
tribution function is predicted. 
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